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Abstract 

A statistical model of quantization based on an exponential distribution of infinitesimal action 



is proposed. Trajectory which does not extremize the action along an infinitesimal short segment 
of path is allowed to occur with a very small probability following an exponential law. Planck 



constant is argued to give the average deviation from the infinitesimal stationary action. 
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I. MOTIVATION 



As early as 1926 Madelung has shown that if one writes the complex-valued wave func- 
tion into polar form, then the Schrodinger equation can be decomposed into a modified 
Hamilton- Jacobi equation and a continuity equation for a probability flow jl|. This obser- 
vation has inspired many efforts to develop statistical models that lead to the derivation 



of the Schrodinger equation 



-|l8|. It also provides the basis for the development of inter- 



pretation of quantum theory in term of traditiona 



include the hydrodynamics interpretation 



QBE 



classical statistical mechanics. These 



25j and Bohmian mechanics 26|-|3l|. Let 



us note however that the idea that the velocity of the particle is determined by the phase of 
a quantum wave in Bohmian mechanics is firstly developed by de Broglie as early as 1923. 



32] 



See for example the interesting historical description in Ref. 

Further, despite the prominent role of Planck constant in the canonical quantiza- 
tion through replacement of c-number (classical number) with q-number (quantum num- 
ber/Hermitian operators), its physical and dynamical origin, after more than one century 
since its first identification which signifies the birth of quantum theory, is still not clear. 
This question is of course beyond the standard formalism of quantum mechanics and can 
only be discussed if quantization is shown to arise effectively from a deeper theory. Hence, 
the elucidation of the physical origin of canonical quantization and thus of Planck constant 
might be indispensable for the searches of violations of quantum mechanics 33H44|. and 



may also lead to useful physical insights in the attempts to search for new physics within 
Planck scale, in which gravitational effect is no more ignorable. 
On the other hand, in the previous work 



45] , we have developed a quantization method for 
systems of spin-less particles based on replacement of c-number by c-number parameterized 
by an unbiased random variable A to modify the classical Hamilton- Jacobi and continuity 
equations. The resulting modified Hamilton- Jacobi and continuity equations can then be 
rewritten into the Schrodinger equation when the distribution of A takes the form 

P(\) = ±6(\-h) + ±6(\ + K). (1) 

Hence, A is an unbiased binary random variable which can take values only ±h. We then 
read-off a unique quantum Hamiltonian from the Schrodinger equation. Unlike canonical 
quantization, the method does not suffer from the problem of operator ordering ambiguity 
and the quantization processes (replacement) has direct and explicit interpretation as a 



specific modification of classical dynamics of ensemble of trajectories (in configuration space) 
parameterized by an unbiased random variable A. 

In Ref . 0] , we have attempted to provide an argumentation that the rules of replacement 
postulated in Ref. ^5] can be derived from a Hamilton- Jacobi theory with a specific random 
constraint uniquely determined by the classical Lagrangian. In the present pap er, we shall 
give an alternative derivation of the rules of replacement postulated in Ref. 45[ from a dis- 
tinct statistical model based on the assumption of exponential distribution of infinitesimal 
action. Trajectory which does not extremize the action along an infinitesimal short seg- 
ment of path is allowed to occur with a very small probability following an exponential law 
characterized by Planck constant. The Schrodinger equation and canonical commutation 
relation are shown as the implication of the model rather than postulated. Moreover, the 
configuration of the system evolves continuously in time and its effective velocity is related 
to the wave function in formally the same fashion as in Bohmian mechanics 
unlike the latter, the wave function is not physically real. 
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3l|. Yet, 



II. EXPONENTIAL DISTRIBUTION OF INFINITESIMAL ACTION 

Let us assume that in microscopic time scale the dynamics of the system is effectively 
stochastic. Further let us assume that the velocity of the system depends on a random 
variable A, q = dq/dt = q(q,X;t). Let us then consider an action along a segment of path 
connecting two infinitesimally close (configuration) spacetime points as: 

rq(t+dt) 

X(A)= / L(q,q(q,X;t))dt, (2) 

Jq(t) 

where L is the classical Lagrangian which now depends on A through q(q,X;t). Let us as- 
sume that A is fixed along the short segment of path. This allows us to apply the principle 
of stationary action in the usual fashion which prescribes that the admissible classical tra- 
jectory connecting the two points is the one which extremizes the action. Solving <5X(A) = 
with variation that vanishes at the end points, one obtains the Euler-Lagrange equation 
(d/dt)(dL/dq) — dL/dq = 0. Let us assume for simplicity that the Lagrangian is not singu- 
lar, det{<9 2 L/dqidqj} ^ 0. Then the Euler-Lagrange equation can be directly rewritten into 
the Hamilton equation 

q{q,p(qA;t);t) = p(q,X;t) = ( 3 ) 
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where p(q, q(q, A; t);t) = dL/dq is the conjugate canonical momentum and H_(q,p) = pq — L 

is the classical Hamiltonian which now depends on A through p(q, A; t). 

Let us proceed to go to the Hamilton- Jacobi formalism through a canonical transforma- 

□ 

tion so that the new classical Hamiltonian is vanishing. One then gets |47J 

-H{q,p) = d t S, (4) 

where S_(q, A; t) is the generating function of the canonical transformation satisfying 

P = d q S. (5) 

Inserting Eq. ([5]) into the left hand side of Eq. (T4]) one obtains the Hamilton- Jacobi equation 
d t S + H((l,d q S) = 0. 

Equations @ and ([5]) are equivalent to dS_ = d t Sdt + d q S_ ■ dq = —Hdt + p ■ dq = Ldt so 
that dS_ is the infinitesimal action along the classical (stationary) path. Now let us define a 
new function S(q, A; t) so that the difference along an infinitesimal stationary path is given 
by dS = dS_ = Ldt = —Hdt + p ■ dq. Using S, the classical dynamics of the system along 
an infinitesimal segment of path with a fixed value of A can then be written in probabilistic 
form as: the probability that an infinitesimal segment of path occurs with dS is given by 

P s (dS\dS) ~ 5{dS - dS) = 5(dS -(p-dq- Hdt)), (6) 

where 5 now refers to the delta function. 

Written in the above form, one can then statistically modify the classical dynamics by 
allowing dS to fluctuate around dS_ = p ■ dq — Hdt. Let us then assume that there is an 
infinitesimal (microscopic) time scale dt = tq so that the system can take a non-classical 
(non-stationary) segment of path with a value of dS, whose probability of occurrence is 
determined by its deviation from dS_ according to an exponential law as follows: 

P s (dS\dS) ~ e -Ws-ds) e -e{s) d t^ (7) 

where 9{S) is a function of S evaluated at the initial point of the segment of the trajectory, 
whose form, as will be discussed below, is determined uniquely by the classical Hamiltonian. 
Evidently, A can not be vanishing and has to have the dimension of action. We assume that 
dS > dS_ for A > and conversely, when A < 0, one assumes dS < dS_. Classical dynamics 
is thus regained in the regime when \dS_/X\ > 1 or formally when |A| <C 1. In this case, 
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Ps(dS\dS_) approaches a delta function centered at dS = dS_ of Eq. (jSJ). |A| thus gives the 
average deviation of dS from dS_. 

Let us discuss how the assumption put in Eq. (J7J) modifies pair of Eqs. (BJ and (JHJ). Let 
us first denote the joint-probability density that the configuration of the system is q with 
the value of the random variable A at time t by fl(q, A; t). The marginal probability densities 
are then given by 

p{q;t) = J d\tt k P(A) = J dqVt, (8) 

where we have assumed that the probability density of A is stationary (independent of time). 

Given a fixed value of A, let us consider two infinitesimally close spacetime points (q; t) 
and (q + dq; t + dt). Let us assume that for this value of A, the two points are connected to 
each other by a segment of trajectory with dS = dtSdt + d q S ■ dq. Then, for a fixed value of 
A, according to the conventional probability theory, the probability density that the system 
initially at (q;t) traces the segment of trajectory and end up at (q + dq;t + dt), denoted 
below as Cl((q + dq,X;t + dt) | \(q, A; t)), is equal to the probability that the configuration 
of the system is q at time t, Q(q,X;t), multiplied by the probability of occurrence of the 
segment of trajectory which is given by Eq. (JTj). One thus has 



tt((q + dq,\;t + dt)\\(q,\;tj^ ~ £l(q,\;t) 
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xe -f(dS-dS) e -0(S)dt_ ^ 

Expanding the exponential on the right hand side up to the first order one gets £l((q + 
dq,X;t + dt)\\(q,X;t)) w [l - (2/X)(dS - dS) - 6(S)dt]tt(q, X;t). This can be rewritten as 

dQ(q,X;t) = - \-(dS - dS) + 6{S)dt\tt{q, X;t), (10) 
LA J 

where dVt{q, A; t) = Q((q + dq, X;t + dt)\\ (q, A; t)) — fi(g, A; t) is the change of the probability 
density Q(q, A; t) due to the transport along the segment of trajectory. 

Taking the limit S — >■ S_, Eq. (fTUj) reduces into dfl(q, A; t) = — Q(q, A; t)9(S_)dt. Dividing 
both sides by dt and taking the limit dt — > 0, one obtains 

^ + e(S)n = d t n + q ■ d q n + e(s)n = o. (11) 

To have a smooth correspondence with classical dynamics of ensemble, the above equation 
which describes the ensemble of classical trajectories, has to be equivalent to the classical 
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continuity equation given by 



6(S) = d q -q(S) = d q -(^ ). (13) 

P=O q S_/ 



d t n + d q ■ (q(S)tt) = o t n + q (S) ■ d q n + d q ■ q (S)n = o, (12) 

where the functional form of q with respect to S_ is given by substituting Eq. ([5]) into the 
left equation of ()3]). Comparing Eqs. (fTT|) and ffl2]) . 9(S_) thus has to be identified as the 
divergence of the classical velocity field 

dH 
V dp 

Accordingly, it is sufficient (while not necessary) to assume that the functional form of 6(S) 
in Eq. ([7]) is given by replacing S_ in Eq. ( |T3l) with S 

6{S) = d q -( d -= ). (14) 
V dp p=d q sJ 

Let us go back to Eq. (jTUI) . Writing dfl and dS as dF = d t Fdt + d q F ■ dq, recalling 
dS_ = d t Sdt + d g S_ ■ dq = —Hdt + p ■ dq, and comparing term by term one finally obtains 

P(Q, O) = d q S(q, A; t) + - 
-H(Q,P(Q,0)) = d t S(q,X;t) + ~ + ^O(S). (15) 



This is just the rules postulated in Ref. 



451 ] where 



ormal "replacement" there is shown here 



as physical "substitution" (see Eq. (5) of Ref. j45|). It is evident that as expected, in the 
formal limit A — > 0, one regains ([5]) and (jlj) respectively. The above pair of equations can 
thus be regarded as the generalization of Hamilton- Jacobi theory. 

To see how Eq. ( |T5l) modifies the classical dynamics of ensemble of trajectories, one thus 
needs to combine it with Eq. f|T2|) . Inserting the upper equation of f|T5|) into the left hand side 
of the lower equation, one obtains a modified Hamilton- Jacobi equation. On the other hand, 
inserting the upper equation of (TT5]) into Eq. f|T2l) one gets a modified continuity equation. 
This is already shown in Ref. (451 to reproduce the results of canonical quantization for a 
wide class of classical systems of point-like particles with no spin if Q(q, A; t) = p(q, |A|; t)P(X) 
with -P(A) is assumed to have the form given by Eq. (CQ). 

Next, inserting the first equation of (fT5l) into the left equation of ([3]), one gets 

\ d 

q = q ( q ,n t) = q(q, d q S(q, A; t) + -^-;t). (16) 
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For example, if the classical Hamiltonian takes the form H(q,p) = p 2 /(2m) + V(q) one 
has g(A) = p/m = d q S/m + (X/2mQ)d q Q. We shall show in the next section that for 
the general case when the classical Hamiltonian is at most quadratic in momentum, 9(S) 
given in Eq. (114p can be interpreted as the divergence of an effective velocity defined 
by averaging g(A) over pair of opposite signs of A. The exponential term e~ 9 ^ dt in Eq. 
(J7J) thus describes whether the segment of trajectory effectively repels {6 > 0) or attracts 
(6 < 0) the nearby trajectories. Intuitively, the probability of occurrence of a segment of 
trajectory that effectively repels (attracts) the nearby trajectories is lower (higher). The 
probability of occurrence of a segment of classical trajectory with dS = dS_ is thus given by 
P s (dS\dS) ~ e - e & dt . 



III. A WORKED EXAMPLE: HAMILTONIAN WITH AT MOST QUADRATIC 
IN MOMENTUM 

For a concrete example, let us apply the above general formalism to a single particle 
subjected to external potentials so that the classical Hamiltonian takes the following form: 

H(Q,P) = ^(p, - Ai)(2. - Aj) + V, (17) 

where Ai, i = x,y,z and V(q) are vector and scalar potentials respectively, g % ^{q) may 
depend on the configuration, and summation over repeated indices are assumed. Inserting 
into the left equation of ([3]) one gets 

q\q^t) = ^ = g l Hp_ r A ] ). (18) 



Substituting this into Eq. (1121 . one has 

d t n + d«(tfi(p.-A j ))ri) =0. (19) 
On the other hand, from Eq. fflBl . 9(S) of Eq. fTH|) is thus given by 

6{S) = d 9l gV{d 9i S-Aj). (20) 
Using the above equation, Eq. ffl5l) becomes 

xd a n 



2 = 9 « S+ 27P 
-H(q,p) = a,S + ~ + \d, t 9"(d„S - A,). (21) 



Now let us combine Eq. ( IT9|) and pair of equations in (I2T1) . Substituting the first equation 
of Eq. d2U> into Eq. (USD, one has 

d t n + d qi (g ij (d q .s- A y )n) + ^*a/%n) = o. (22) 

On the other hand, inserting the upper equation of (12"T|) into Eq. (1171) and plugging into the 
left hand side of the lower equation of ( 12~TT) . one has, after arrangement 

dtS + ~ Ai)(d qj S - A,) + V 

2 ^ i? + 9 ^ i? y 

+ (^(^5-^)0) + =o, (23) 

where we have defined i? = y/Ti and used the identity: (1 / 'A£l 2 )d qi £ld q .Q = (l/2Vt)d qi d q .£l — 
(l/R)d qi d qj R. Inserting Eq. fT22|) . the last line vanishes to give 

dtS + ^(d qi S - A){d qj S - A,) + V 

X 2 / -da d a R .d a .R\ , 

-TK 9 ^r~ + d « 9 r) =0 - (24) 

The above equation reduces into the Hamilton- Jacobi-Madelung equation when g l i is inde- 
pendent of q and X = ±h. We have thus the pair of coupled equations (122"]) and ( 1241) which 
are parameterized by A. 

Let us proceed to assume that f2(g, A; t) has the following symmetry: 

n(q,X;t) = n(q,-X;t), (25) 

so that the probability density of A is unbiased P(X) = f dqQ(q, A; t) = P(— A). In this case, 
S(q,X;t) and S(q,—X;t) + S (X), where S (X) is independent of q and t, satisfy the same 
differential equation of (IM1) : namely the last term of Eq. (|24j) is not sensitive to the signs 
of A. Hence, if initially one has S(q, A; 0) = S(q, —A; 0) + So (A), then one will have 

S(q,X;t) = S(q,-X;t) + S (X). (26) 



So(A) is thus an odd function of A, Sq(— A) = —Sq(X) [48(]. The above properties can then 
be used to eliminate the last term of Eq. ( 1221) : taking the case when A is positive add to it 
the case when A is negative and dividing by two, one gets 

d t Q + d qi [g ij {d qj S - Aj)ri) = 0. (27) 



We have thus the pair of coupled equations (124|) and (j27|) which are still parameterized by 
A. 

Next, since A is non- vanishing, one can define the following complex- valued function: 

v|/(g,A;t) = i?exp(^s). (28) 

It differs from the Madelung transformation in which 5* is divided by |A| instead of h. The 
probability density for q is thus given by 

p(q-t) = j d\Q = j ' d\\^\ 2 . (29) 

Equations fT2"4"j) and f|27|) can then be recast into the following modified Schrodinger equation: 

i\\\d t ^ = ~H|A|d 9i - A i )^'(g)(-i|A|5 ? , - A,-)* + VV. (30) 

Here we have assumed that the fluctuations of A in space and time are ignorable as compared 
to that of S. Let us then assume that Q is factorizable as fl(q, A; t) = p(q, |A|; t)P(X), where 
P(A) takes the form given by Eq. (JT]). In this case, Eq. (1301 reduces into the celebrated 
Schrodinger equation 

m t ^ Q {q-t) = HV Q (q;t), 
* Q (q; t) = y/p(q,h;t)e* s ^ & S Q ( g] t) = S(q, ±H) t), (31) 

where the quantum Hamiltonian H is given by 

H = ^ - A i )g i ^q)(p j - A,) + V, with ^ = -ihd qi . (32) 

From Eq. (13T1) we know that the Born's statistical interpretation of wave function is valid by 
construction for all time, p(q;t) = \^Q(q; t)\ 2 . In this context, Eq. ( 1291) should be regarded 
as a generalization of Born's rule. One can also see that unlike canonical quantization, the 
present statistical model of quantization selects a unique operator ordering in which g^{q) 
is sandwiched by p — A. 

The Schrodinger equation is thus a subset of the present model corresponding to a discrete 



unbiased random variable A = ±h. It is then interesting to mention Ref. 49] that the master 
equation of a particle moving with a fixed velocity, subjected to a random complete reverse 
of direction following a Poisson distribution, can be written into a Dirac equation (in the 
same way that the Schrodinger equation is connected to the dynamics of Brownian motion) 



through analytic continuation. Further, while A = ±h is a binary random variable, it might 
be a function of a set of continuous "hidden" variables £ = {£1,^2, £3, ■ • • }: A = /(£). For 



example, one may assume that A = \JQ + £| + £f = ±ft. Hence, £ = {£1,^2, £3} lies on the 
surface of three-dimensional ball of radius h. Let us divide the ball into two with equal area 
and attribute ±h to each division. If £ moves on the surface sufficiently chaotically, then 
one will obtain A = ±H with equal probability. 

Since Q and S are multi- valued functions of q due to their dependence on A, then ^(q, A; t) 
is also a multivalued function of q. Let us consider the case when A can only take binary 
values A = ±h as assumed in Eq. ([Q). Equations ( 1251) and ( 1261) then become 

Q(q, H) = Q(q, -K) & S(q, K) = S{q, -K) + S {h). (33) 

Hence, the probability density becomes single-valued while the phase is still multi-valued in 
accord with quantum theory. In contrast to quantum theory, the wave function is however 
multi- valued satisfying 

y Q (q,H) = y Q (q,-h)e^ h \ (34) 
Assuming that the wave function is single-valued ^Q(q, K) = ^Q{q, —h) one then has 

S = nh, n = 0,±l,±2,..., (35) 

where h = 27rH is the Planck constant. Hence, the set of single-valued quantum wave 
functions is a subset of the wave function of the present model when one initially picks a 
phase satisfying S(q, h; 0) = S(q, —h; 0) + nh, n = 0, ±1, .... 



Wallstrom then argued in Ref. [50] that the Schrodinger equation with multi- valued wave 
functions does not give discrete quantum number, say discrete values of angular momentum. 
As discussed in Ref. |45|, in our model, a system can indeed possess any continuum values 
of (effective) angular momentum, energy etc. As shown there, however, a measurement of 
angular momentum can only give discrete possible values, that is one of the eigenvalue of 
quantum angular momentum operator as in standard quantum theory. Namely, quantum 
discreteness is a feature of measurement. 

Further, inserting the upper equation of (12T]) into Eq. (TT%|) one has 

A d ^1 

q\\) = 9 %] {d q3 S~A J ) + -g^^ r . (36) 
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One can then define an effective velocity as 



m = m+ f~ X) =3"(d, 1 S-A 1 ), (37) 



where the second equality is obtained by the virtue of Eqs. (I25j) and (126|) . Putting A = ±H, 
it reduces into 

v\h)= g v(d q] S Q -A j )=v\-h). (38) 

This is just the actual velocity of the particle in Bohmian mechanics, satisfying the continuity 
equation 

d t p + d q ■ (pv) = 0. (39) 

Hence, we have an effective picture that the particle moves as if it is guided by the wave 
function a la Bohmian mechanics (de Broglie-Bohm pilot- wave theory). However, unlike the 
latter, the wave function in the present statistical model is not physically real. It is just an 
artificial mathematical tool to describe the dynamics and statistics of ensemble of copies of 
the system. Let us also mention that the effective velocity field v % of Eq. fl38l) is equal to the 



"naively observable velocity field" reported in Ref. 51] obtained through weak measurement 



52j within the standard formalism of quantum mechanics. 



IV. CONCLUSION AND DISCUSSION 

We have thus argued that quantization may be understood as the implication that the 
principle of stationary action is valid only approximately, and the distribution of the devia- 
tion from the stationary action along infinitesimal segment of path for a certain microscopical 
time scale tq follows the exponential law given in Eq. (JTj). Moreover, in this statistical model, 



the absolute value of the random variable | A| = h postulated in Ref. [45] is interpreted as the 
average deviation from the infinitesimal stationary action. The origin of Bell-non-locality 
in many particles system may then be traced back to the fact that dS is evaluated along 
a segment of trajectory in configuration space rather than in ordinary space. Namely, a 
probability density is attributed to the infinitesimal action along a segment of trajectory in 
configuration space, rather than to the trajectory of each particle in ordinary space. 

While we have argued that the reduced Planck constant h is the average deviation from 
the infinitesimal stationary action in a certain microscopical time scale tq, the model does 
not offer physical mechanism that determines the value of the Planck constant. An attempt 
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;o provide a cosmic origin of the value of Planck constant is proposed by Calogero in Ref. 
53l | . Let us also mention the work by de la Pena and Cetto reported in Ref. 54| which 
is partially motivated by the work of Calogero attempting to explain the value of Planck 
constant within the stochastic electrodynamics 5^|. 

Nor the present statistical model explains why the deviation from the infinitesimal sta- 
tionary action follows an exponential law given by Eq. ([7]) rather than other distributions. 
To discuss this issue in the context of our model, it is then necessary to investigate the 
physics beyond the quantum mechanical time-scale tq. In particular, it is interesting to ask 
if one can devise a deterministic dynamical system which is valid in time scale less than 
tq, leading effectively in the quantum mechanical time scale of tq to a stochastic behavior 
governed by Eq. ([7]) rather than other distributions. A proposal that quantum mechanics 
is emergent from a deterministic model for example is given in Ref. |56|. To this end, it is 
also interesting to study the relation between the statistical model presented in this paper 
and that developed in Ref. 46[ based on Hamilton- Jacobi theory with a specific random 
constraint determined uniquely by the Lagrangian. 

The present statistical model suggests several possible generalizations of the Schrodinger 
equation of Eq. fl3T|) . The first immediate possible generalization is given by Eq. f )30|) 
which is valid even when the distribution of A deviates from Eq. ([1]) but is still unbiased 
-P(A) = P(— A). For example, one assumes that A fluctuates around ±h with a very small 
width. Such type of generalization is discussed in Ref. 5_3j and is shown to lead to testable 
possible corrections to the statistical prediction of quantum mechanics. The other type 
of generalization is to choose different form of 0(S) which still reduces into Eq. ( 1T31) in 
the limit S — > S_ to guarantee smooth classical correspondence, and gives Eq. (TH|) as 
its first order approximation. Yet another possibility is to take different distribution of 
infinitesimal action which reduces into Eq. in some physically reasonable limit. All the 
above mentioned problems together with the application of the present statistical model to 
foundational problems of quantum mechanics will be the object of future investigation. 
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